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Abstract 

We consider back-reaction by non-supersymmetric D7/anti-D7 probe branes in the 
Kuperstein-Sonnenschein model at finite temperature. Using the smearing technique, we 
obtain an analytical solution for the back-reacted background to leading order in Nf/N c . 
This back-reaction explicitly breaks the conformal invariance and introduces a dimension 6 
operator in the dual field theory which is an irrelevant deformation of the original conformal 
field theory. We further probe this back-reacted background by introducing an additional set 
of probe brane/anti-brane. This additional probe sector undergoes a chiral phase transition 
at finite temperature, which is absent when the back-reaction vanishes. We investigate the 
corresponding phase diagram and the thermodynamics associated with this phase transi- 
tion. We also argue that additional probes do not suffer from any instability caused by the 
back-reaction, which suggests that this system is stable beyond the probe limit. 
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1 Introduction 



The understanding of the gauge-gravity duality or the Anti-de-Sitter/Conformal Field Theory 
(AdS/CFT) correspondence [U El E] (see e.g. jl] for an earlier review) has provided us with a 
large class of strongly coupled large N c gauge theories. The very nature of this duality enables 
us to perform computations in classical general relativity which — when translated to the gauge 
theory language via the AdS/CFT dictionary — corresponds to a computation in the dual field 
theory in the strong coupling limit. This is a remarkably powerful technique to probe aspects of 
strongly coupled gauge theories where conventional perturbative field theory methods are severely 
inadequate. 

Strongly coupled systems are abundant in nature: the physics of the strong interaction de- 
scribed by Quantum Chromodynamics (QCD), as suggested by its name, describes a strongly 
coupled system. Coincidentally we live in an interesting time when the strong-coupling physics of 
QCD has been experimentally explored in the Relativistic Heavy Ion Collider (RHIC) and is being 
currently explored in the Large Hadron Collider (LHC); see e.g. [HJ El El El E] for more information 
about experimental results and findings. These explorations and findings necessitate improved 
theoretical control in addressing strong coupling physics. 

This is where the AdS/CFT correspondence makes heiQ entrance. Currently we do not have 
a gravitational dual to QCD in a precise sense. However, within this class of large N c gauge 
theories, the strongly coupled dynamics of the "quarks" and the "gluons" can be explored and 
understood. Such theories are typically supersymmetric and conformal with adjoint matter fields 
(= "gluons"), and in addition fundamental matter (= "quarks") can be introduced. The hope is 
that by studying these theories, we can learn general lessons about universal features of strongly 
coupled systems that also apply — at least qualitatively — to QCD. 

In this article, our focus will be both on the fundamental flavor sector and its back-reaction 
on the adjoint sector. Typically the fundamental matter fields are introduced by considering Nf 
number of flavor branes in the background of N c color branes. The stack of color branes gives 
rise to a gravitational background which can be obtained by solving the supergravity equations of 
motion. In this background the flavor branes are introduced in the probe limit, with Nf <C N c , 
such that they do not back- react on the background. This idea was pioneered in (TO]. Many 
interesting physical questions can be explored and understood in the flavor sector within this 
probe limit and within various models, e.g. the meson spectrum and the meson melting phase 
transition respectively in [TT] and [T21 [T5| H3] in the D3-D7 model. Another canonical example is 
the physics of the chiral symmetry breaking within the flavor sector realized in e.g. the so called 

1 We apologize if this gender assignment appears presumptuous. It seems the age-old folklore of a conventional 
assignment to the explorer and the explored has gotten the better of us. 



1 



Sakai-Sugimoto model [15j [16] and its finite temperature version in [T7J [IH]J^] and more recent 
models of chiral symmetry breaking in [20l |2T] F] 




All these analyses are performed within the probe limit where it is assumed that the number 
of flavors is very small compared to the number of colors. In QCD, however, this is not the case. 
First we have N c = 3 and thus large N c gauge theories will not work as a good quantitative 
approximation and we may only hope for robust qualitative features. Second, we have Nf ~ iV c ; 
hence the physics that we observe in the strict probe limit may also be of limiting use. One may 
wonder how much the qualitative physics can really depend on the precise values of Nf and iV c 
and their relative strength. It certainly depends on the precise physical question, however we 
can identify certain physical properties that crucially depend on these values. One such example 
is the QCD beta function: it can be shown that QCD is an asymptotically free theory as long 
as Nf < yiV c , which holds true in nature. However, the physics will be completely different if 
the above inequality is not satisfied. Another crucial property is chiral symmetry breaking where 
the physics of the chiral transition depends on N c and Nf leading to a phase structure which is 
sensitive to these numbers (at least on their ballpark ranges) QThus we can safely conclude that 
going beyond the probe limit, as far as introducing flavor degrees of freedom is concerned, is a 
rather important topic to explore]^] 

The general procedure to consider back-reaction consists of the following steps: We start 
with a background obtained by solving the supergravity equations of motion, e.g. an AdSs x 
M 5 background, where M 5 is some Sasaki-Einstein manifold. This gravitational background is 
obtained from the closed string sector. The dual field theory is some M = 1 superconformal 
quiver gauge theory with a certain global symmetry group which is identified with the isometry 
group of the internal manifold. To introduce flavors, we need to add open string degrees of freedom. 
The dynamics of the open string degrees of freedom is described by the probe D-brane action of 
appropriate dimension and is given by the so-called Dirac-Born-Infeld (DBI) action (supplemented 
by the Wess-Zumino term whenever necessary). To consider the back-reaction we need to consider 
the DBI action as the matter source that deforms the original background. 

Work along this direction has been undertaken by several authors over the years beginning 
with e.g. [2U I25J [213 [2H [2E1 [2H] and nicely summarized in [3D]. These references have analyzed 
the back-reaction by supersymmetric probes at zero temperature. Supersymmetry makes the 
problem more tractable since one does not need to solve the equations of motion, but instead it is 

2 It is argued in |19| that the chiral symmetry restoration mechanism proposed in |171 118) is problematic and a 
new mechanism is proposed. We thank Takeshi Morita for pointing this out to us. 

3 For a comparative study of the D3-D7 model and the Sakai-Sugimoto model in the presence of various external 
parameters, see e.g. [22] . 

4 See e.g. [33J for some lattice results. 

5 This limit is known as the Veneziano limit. 
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sufficient to solve first order BPS equations. This technical simplification comes at the cost that the 
resulting theory — even after the inclusion of the flavor degrees of freedom — is supersymmetric 
and does not realize the spontaneous breaking of chiral symmetry as one would expect in QCDjj 




Nonetheless, numerous interesting features have been explored within this type of models at zero 
and finite temperature and more recently including other parameters, such as a chemical potential 
or an external electromagnetic field in e.g. [32], [331 EH EH |36l [37] . 

Here we will focus on a model recently discussed in [20] where D7 and anti-D7 branes have 
been introduced in the AdSs x T > background. The brane-anti-brane pair wraps a three cycle in 
the internal manifold T 1 ' 1 = S 2 x S 3 and extends along the rest of the cone. At zero temperature, 
the brane-anti-brane pair joins in the IR realizing a spontaneous breaking of the chiral symmetry: 
U(Nf) L x U(Nf)n — > U(Nf) diSbg . The asymptotic angle separation between the brane and the 
anti-brane — previously denoted by A^, but will henceforth be denoted by A0* for reasons that 
will become apparent in the main text — takes a constant value and corresponds to the coupling 
of the corresponding operator in the flavor sector. In [38J we have discussed in details the phase 
structure of this model in the presence of a non-zero temperature and an external electro-magnetic 
field. 

Before going further, let us comment on the relevance of the Kuperstein-Sonnenschein model. 
First, even at zero temperature, the probes break supersymmetry completely and realize a spon- 
taneous breaking of chiral symmetry. Furthermore, unlike the Sakai-Sugimoto model, the dual 
field theory is an honest (3 + l)-dimensional one. Thus, this model potentially realizes certain 
properties of QCD better than the Sakai-Sugimoto model. In this article, we will consider the 
back-reaction by the flavor sector in this model at finite temperature. 

The rationale behind considering back-reaction in the finite temperature case is three-fold: 
first note that since the probes break supersymmetry completely, no simplification such as solving 
first order BPS equations will occur at zero temperature. Second, the profile at zero temperature 
is described by a non-trivial function 0(r) (where <fi is the azimuthal angle of the S 2 where the 
probes are extended), which makes analyzing the back-reaction problem more involved]^] Third, 
for our current purposes — which is to analyze the phase structure once the back-reaction is taken 
into account — the relevant background should have a non-zero temperature. 

The astute reader might raise the following question: how is this system any different from what 
has been obtained in e.g. [32] which also deals with a finite temperature back-reacted solution? 
This is indeed a very relevant question. For now we will offer a brief answer to this question 

6 Note that the issue of the back-reaction in the Sakai-Sugimoto model, where both supersymmetry is completely 
broken and we have a geometric realization of spontaneous chiral symmetry breaking, has been addressed in |31j . 

7 Recall that when the background is AdS-SchwarzschildxT 1,1 , the favorable embeddings are described by <fi = 
const as shown in |38) . 
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through some technicalities involved in the computations and elaborate on these throughout the 
rest of the paper. First recall that the precise embedding as considered in e.g. [25] and the ones 
that we consider here a la [20] are different in the details. These details turn out to be important 
for the physics that we eventually extract from such back-reacted systems. Second, in our system, 
we have both D7 and anti-D7 branes as opposed to having just branes in the systems considered so 
far. The D7-branes are codimension two objects and therefore they act as two dimensional delta 
function sources. The resulting equations of motion become coupled PDEs which pose a challenge. 
This difficulty can be circumvented by using the so-called "smearing" technique which means we 
introduce many such probes and distribute them in the transverse directions homogeneously. If we 
have an equal number of branes and anti-branes smeared at each point, it will not source any Cs 
potential (the Hodge dual to the axion field). This is in stark contrast to the framework considered 
in e.g. [32] where there are no anti-branes and therefore a non-trivial Cg potential is sourced. 

Complying with the conventional wisdom, we find that the flavor back-reaction is understood 
as irrelevant deformations of the original CFT. In our case, these deformations consist of the 
departure from conformality and the insertion of a dimension 6 operator. The former originates 
from a running dilaton and the later from the deformation of both the S 2 x S 2 base and the U(l) 
fiber direction of the original T 1,1 -manifold. We will not make attempts to precisely identify this 
dimension 6 operator, however we will extract a lot of information about the "phase structure" of 
this theory in a similar spirit as in [38]. 

The breaking of the conformal invariance due to flavor back-reaction has a very important 
physical consequence. To appreciate this, let us recall that in the original Klebanov-Witten back- 
ground the flavor sector did not undergo a chiral phase transition at finite temperature simply 
because there was no other scale available [38J. The back- reaction of the flavors introduces this 
dimension 6 operator, equivalently a dimension —2 coupling denoted by m~ 2 , providing another 
natural scale in the system. We can now probe this back-reacted background with an additional 
set of N'x D7/anti-D7 branes and focus on the physics this extra set of flavor sector experiences. 
The interplay of the background temperature and the dimension —2 coupling conspires to conjure 
up a non-trivial phase structure. We have studied this phase structure in the A0* vs m 2 plane, 
where m 2 is defined as 



Here T is the background temperature, L is the radius of the original AdS space and A is the 't 
Hooft coupling. 

Another important and interesting question regarding such non-supersymmetric probing is the 
issue of stability: it is not entirely clear how stability might be affected once the back-reaction 
by the non-supersymmetric probes are taken into account. It was demonstrated in [39J that the 
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non-supersymmetric probes are stable in the Kuperstein-Sonnenschein model. Here we will not 
undertake the task of analyzing the fluctuation modes of the back-reacted background. We will, 
however, analyze the stability of the additional set of probes that we introduce in this back-reacted 
background. We will argue that this additional probe sector does not experience any unbalanced 
force due to the absence of supersymmetry and the fluctuations of these additional probes are 
not likely to experience any instability caused by the back-reaction itself. Thus this is also very 
suggestive of the fact that perhaps stability is not a matter of concern in this model in which 
supersymmetry plays no role. We will comment more on this issue in the concluding part of this 
work. 

This paper is organized as follows: we start with a brief review of the Kuperstein-Sonnenschein 
model at finite temperature in section 2. In section 3, we discuss the issue of taking back-reaction 
into account, describe the analytical solution treating the back-reaction perturbatively and offer 
comments on the dual field theory. In the following section, we discuss various thermodynamic 
properties of this back-reacted background. Then, in section 5, we introduce an additional set of 
probe brane/anti-brane system in this back-reacted background and study the phase structure in 
the probe sector. In section 6, we argue that this probe sector does not suffer from any instability 
after the inclusion of the back-reaction. Finally we conclude in section 7 with future directions. 
Various technical details have been relegated to four appendices: appendix A contains a discussion 
of the smearing form, appendix B contains the details of the equations of motion, in appendix 
C we have provided the detailed steps towards obtaining the leading order back-reacted solution 
and in appendix D we have obtained the action corresponding to the quadratic fluctuations of the 
additional probe sector in our back-reacted background. 

2 Kuperstein-Sonnenschein probes at finite temperature 

In [20], the authors introduced probe D7/anti-D7 branes in the AdSs x T 1 ' 1 background such that 
the D7/anti-D7 branes wrap a three cycle in the internal manifold T 1 ' 1 = S 2 x S 3 and are extended 
along the rest of the conifold IR 2 x S 2 . The field theory dual to this supergravity background is 
given by an A^ = 1 superconformal quiver gauge theory with a gauge group SU (N c ) x SU (N c ) and 
two bi-fundamental chiral superfields which are usually denoted by and Bi 2 - These chiral 
superfields transform in the (iV c , N c ) and the (N c , N c ) representations of the gauge group. This 
theory also possesses a global SU{2) x SU{2) x U{1)r symmetry which, in the dual gravitational 
background, manifests itself as the isometry group of the internal manifold T 1 ' . 

Introducing probe D7/anti-D7 branes amounts to introducing matter fields in the fundamental 
representation: specifically it corresponds to adding left-handed/right-handed Weyl fermions in 
the dual gauge theory. This gives rise to a global symmetry group of U(Nf) L x U(Nf) R flavor 
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symmetry, where Nf is the number of flavors that we have introduced in this background. 

At zero temperature, the probe D7 and anti-D7 branes have no choice but to dynamically join 
in the IR and realize spontaneous breaking of chiral symmetry: U (Nf) l x U (Nf) R — > t/(^) diag . In 
the presence of a black hole horizon in the bulk (and therefore finite temperature in the dual field 
theory), there are two types of profile: one where the brane-anti-brane pair joins at some IR radial 
position; and the other where they separately fall into the black hole. The former corresponds 
to a chiral symmetry broken phase and the later to the chiral symmetry restored phase. This 
phenomenon is pictorially represented in fig. [TJ 



Figure 1: A schematic representation of the shapes of the probe branes. Fig. (a) The red dot 
denotes the conifold singularity located at r = 0. The probe branes join at r = ro realizing the 
spontaneous breaking of chiral symmetry. Fig. (b): Two types of embedding are possible at finite 
temperature and the only available phase is the chiral symmetry restored phase. 

Although we have two different types of embeddings, there is no possibility of a phase transition 
between them. The reason is simple: since our background is conformal, temperature is the only 
scale. Thus there is no meaning of having a "high" or a "low" temperature phase, there will be an 
unique phase at any temperature. As argued and explicitly shown in [38j, introducing a non-zero 
temperature immediately favors the chiral symmetry restored phase. 

Let us be more specific. The background at finite temperature is given by AdS 5 - 
SchwarzschildxT 1,1 j^] The time direction is Euclideanized and compacitified on a circle. The 

8 In this case the dilaton is a constant, therefore the string frame metric and the Einstein frame metric are the 
same up to a factor of the string coupling g s . In keeping with the notations used in the later parts of the paper, 
we will define g s := e** . 




(a) 



(b) 
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temperature is then simply given by the inverse period of this circle. In Lorentzian signature this 
background is given by 

ds 2 = — (-b(r)dt 2 + dx 2 ) + + R 2 ds 2 T ^ , b(r) = 1 - -J , (2.2) 

L r j[r) r 



o 



(e 1 ) 2 +(e 2 ) 2 + (g 1 ) 2 +(e~ 2 ) 2 ]+i(e 3 + e 3 ) 2 . (2.3) 



Here t is the time, x represents the spatial 3-directions, r# is the location of the event-horizon and 
L is the radius of curvature of the AdS-space. After Euclideanization, the temperature is given 
by T = r H /(nL 2 ). Also e i and e* are the left-invariant one- forms of two independent SU (2)'s. 
There is an £7(1) symmetry which rotates e 3 to e 3 . Finally there is a Z2 symmetry which acts as 
Z 2 : e 1 ' 2 -H- e 1,2 . Thus the background has an SU(2) x SU(2) x U(l) X Z 2 symmetry, where the 
£7(1) is the R-symmetry. 

We choose the vielbeins as follows [4 lj 

ip ib ib ip 

e 1 = sin — sin 9\d<p\ + cos ~^d9\ , e 2 = — cos — sin Q\d§\ + sin ~^d9\ , (2.4) 

ip ip ip ip 

e 1 = sin — sin 9 2 d(p2 + cos — d9 2 , e 2 = — cos — sin 9 2 d(p 2 + sin — <i6* 2 , (2.5) 

Zi Zi £i Zi 

e 3 + e 3 = + cos #i#i + cos #2^2 • (2.6) 

Here and {#2,02} parametrize the two if> 2 s. The ranges for the angular directions are 

given by: ip G [0,47r), <pi G [0, 2tt) and 9i G [0,7r). With this parametrization, the metric on the 
T 1 ' 1 becomes 

ds 2 T x,x = - (c?0 2 + sin 2 6W0 2 ) + - ( + ^ cos Oidfa J 

= d4 E + (# + A KE ) 2 , (2.7) 

where in the second line the metric on T 1,1 has been written as the U(l) fibration over a Kahler- 
Einstein manifold S 2 x S 2 and Ake yields the Kahler form via J = gL4ke/2. 

The probe branes are aligned such that they span the Minkowski- and {ip, 6^ 2, (^i^j-directions 
and their profile is represented by {#2,1 {f), 02, 1 (?")}• The DBI action in this case is given by 

Sdbi = -t 7 J d 8 ^detP[G] = -M J dtdrr 3 (l + jb(r) (9' 2 2 + sin 2 9 l>2 /2 2 )) ' , (2.8) 

where T7 is the tension of the probe brane and Af = (t7Vr3 87t 2 )/9. This action leads to equations 
of motion which admit two different classes of solutions represented by 

9 l>2 = vr/2 , li2 = 0(r) , (2.9) 
0i, 2 = vr/2 , 12 = const . (2.10) 
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The profiles described by a non-constant 0(r) are called the U-shaped embeddings and they 
correspond to the chiral symmetry broken phase; the other ones are called parallel-shaped and 
correspond to the chiral symmetry restored phase. These two classes of profiles are schematically 
shown in fig. [T] It was shown in [3E] that energetics always favors the chiral symmetry restored 
phase in the presence of a non-zero temperature. In the subsequent sections we will analyze the 
back-reaction on the supergravity background by these parallel-shaped brane/anti-brane profiles 
and obtain the back-reacted background at finite temperature. 

3 The back-reacted background 

From here on we will discuss the effect of back-reaction by the flavor branes at finite temperature. 
The probes in our system are co-dimension two objects and thus they give rise to delta function 
sources in the plane transverse to the seven branes. Thus, the general problem of considering back- 
reaction by the flavor branes will lead to a set of coupled partial differential equations which is 
difficult to solve. To simplify this problem we can make use of the so-called "smearing" technique 
as demonstrated in e.g. [27]. 

The smearing procedure can be understood as uniformly distributing the large number of 
probes along the transverse directions such that the global symmetry of the original background is 
recovered after taking the back-reaction into account. Following the discussion in [27]) we briefly 
elaborate on this idea. The smearing process for Nf D7-branes can be operationally understood 
as the following modification to the Wess-Zumino contribution 



where f2 is a two- form known as the smearing formjj .Ms denotes the D7-brane world- volume and 
Aiio denotes the full ten-dimensional manifold. This two-form determines the distribution of the 
RR-charges sourced by the D7-branes along the transverse directions. Clearly, if we have an equal 
number of uniformly smeared D7 and anti-D7 branes the Wess-Zumino contributions will cancel 
each other which implies that there is no net RR-charge at each point along the directions where 
the smearing procedure has been implemented]^] 

Moreover, due to the smearing procedure the DBI-piece of the probe action — in Einstein 

9 We will find this smearing form in our case is simply given by a constant. Further details on this computation 
have been relegated in appendix A. 

10 Thus the back-reaction in our system will not source any C$ (the dual to the axion field) potential and 
consequently we will not have any Fi RR-form. This is to be contrasted with the case in e.g. [27] where the 
back-reaction sources a non-trivial axion field. 




(3.11) 
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frame — will take the form 
Sdbi - 



r 7 ! d 8 £e*^/^Gs^-T 7 f d l0 xe*^G^\n\ , 

J Ms J Mia 



(3.12) 



So in view of the above discussion, to consider the back-reaction we will work with the following 
action in Einstein frame 



S 



•SsiTGRA + ^flavour 
1 



2 K io 



d 10 x V-Gio 



R-±d M ®d M ®-- A 



F 5 \ 2 



-G 



(i) 



-G 



(2) 



(3.13) 



(1 2) 

where the superscript on the induced metric G 8 is intended to denote the probes which are 
located at #1.2 = tt/2 and <pi^ = const. Also, kiq is related to the ten-dimensional Newton's 
constant, R is the Ricci scalar, Gio is the background metric, $ is the dilaton and F5 is the self- 
dual five-form field strength. The supergravity part of the action scales as g~ 2 and the flavor part 
of the action scales as g~ l . This is because the supergravity part emerges from the closed string 
sector whereas the flavor part originates from the open string sector. Note that we have defined 
the dilaton field such that $ vanishes asymptotically and the constant $*, which sets the string 
coupling via the relation g* = g s := e**, has been factored outj^j] The quantity r 7 is the tension 
of the brane: r 7 = 5f7 1 (27r) _7 a /4 . The factor of 2 in front of the DBI-action comes from the brane 
and the anti-brane which contribute equally. More details on the action and the corresponding 
equations of motion have been relegated to appendix B. 

Now that we know which action to extremize to find the desired solution, we should start 
with an ansatz for the back-reacted background. The full ten-dimensional metric should be a 
deformation of the AdSs-SchwarzschildxT 1 ' 1 background. More precisely — since the D7-branes 
wrap the entire AdS part — the deformation should appear only within the internal T 1 ' 1 -part. 
Recall that the smearing procedure has recovered the full SU{2) x SU{2) x f/(l) x Z 2 global 
symmetry. Thus the internal manifold is still constraint to possess this isometry and — much like 



11 Here we will use the symbol g* to denote the string coupling in the presence of a non-zero back- reaction in 
keeping with the other notations. 



9 



the case considered in [27] — we make the following ansatz 

ds 2 = h- 1/2 (r) \-b(r)dt 2 + dxidx 1 ] + h l,2 (r){ ^— + V (d9 2 + sin 2 Q { d$) 

V i=l,5 



e a/(r) 

+ — — I # + > ' cosfl, d<Pi\ \ , (3.14) 



F 5 = £;(r)/i(r) 3/4 (V A e" 1 A e x ' 2 A e x " A e r + A e* 1 A e^ 1 A e* 2 A e^ 2 ) , (3.15) 
where h(r), b(r), f(r) and g(r) are unknown functions and the vielbeins are defined as 



br x l% x l 2 dt , e xl = h-^dx 1 , e r = h^b-^dr , (3.16) 

e V = -h}l A e f {dip + cosM0i + cos 2 # 2 ) , (3.17) 
3 

e 01 ' 2 = -=h xl *(*dBx, 2 , e^ 2 = ^h 1/4 e 9 sin9 1>2 d(f) h2 . (3.18) 
V6 V6 

We also have a non-constant dilaton $(r). The quantization condition for the five form flux 

' 1 " N c r 3 , (3.19) 



where T3 is the tension of the D3-branes, gives 



k(r)h(r) 2 e^ r)+ ^ r) = 27irg*NJ 4 s = 4L 4 := ^AZ 4 , (3.20) 
with A = 4irg* s N c . (3.21) 



With the ansatz in (3.14) one can now obtain six independent equations for the five unknown 



functions obtained by extremizing the action in (3.13). These equations are shown and discussed 



in details in appendix B. We will not discuss these details here, instead we will move on to the 
analytical solution that can be obtained for the leading order correction in Nf/N c in the next 
section. 

3.1 Perturbative solution 

We will now exhibit the perturbative solution that we obtained treating the back-reaction as a 
"small" effect. Before discussing this solution explicitly, let us be precise about what we mean by 
this "smallness" . From the action we can explicitly check that the back-reaction part of the action 
(which is given by the DBI-piece) is suppressed by a factor of Nf/g s ~ (XNf)/N c as compared to 
the supergravity action. This is expected since the supergravity action scales as N 2 and the probe 



action in (2.8) scales as (XNfN c ). Thus to obtain a perturbative solution the relevant expansion 



parameter, denoted by e, will be e ~ (\Nf)/N c . 
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To fix the numerical coefficient, we look at the right hand side of e.g. equation (B.24) and 
define 



6 



7T 



9* s N f 



3 
2^ 



XN 



(3.22) 



where we have used the definition of A from equation (3.21). We will discuss the back-reacted 
solution at the leading order in e. 

In fact it turns out that we can find an analytical solution in the leading order in e: in appendix 
C we have discussed the details on how to obtain this solution. Here we will only present the final 
result and discuss its properties. The various functions are obtained to be 



b(r) : 
<3>(r) 
h{r) 

e f(r) 
P 9(r) 



1 



4 

r 



In 



1 + 
1 + e 

1 + e 



24 
1 

48 



+ 4m" 



H 



K 1 



H 



m~ 2 r 2 E I 1 - 



H 



m 



-2 



K 1- 



+ 2m- 2 r 2 E I 1 - 



(3.23) 

(3.24) 

(3.25) 
(3.26) 

(3.27) 



Here r# is the location of the event-horizon, L is the radius of the AdS-space, K and E are elliptic 
functions of the first and the second kind respectively and m~ 2 is an undetermined constant of 
dimension length -2 . Thus we have actually found a one parameter family of solutions parametrized 
by this constant. Also, r* is a UV cut-off. The dilaton runs logarithmically and thus blows up at 
r — > oo, which implies that the string coupling also blows up. Thus we need to define this theory 
with a UV cut-off at r = r* where $ — > 0, 
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Let us now discuss some properties of this solution. Note that the solution for the emblackening 
factor b(r) and the warp factor h(r) do not receive any qualitative correction due to the back- 
reaction. At most we can redefine the AdS-radius L by absorbing the e correction to h(r), which 
will define an effective 't Hooft coupling for us. On the other hand, the functions e* and e 9 do 
receive interesting corrections. At first it might seem that the back-reaction changes the asymptotic 
behavior of these functions, but by the very nature of being a perturbative solution this is not 
the case. We will discuss the regime of validity of this solution in the next section where it will 
become clearer. 



12 Recall that we have already factored out the string coupling from the dilaton and therefore the dilaton has to 
approach zero asymptotically. 
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It is evident from the perturbative solution given in (3.23)-(3.27) that the flavor back-reaction 
deforms the original AdSs x T 1 ' 1 -space and this deformation shows up entirely in the internal 
manifold T 1,1 . Such deformations are understood via the GKPW formula^ [3] by looking at the 
asymptotic behavior of the fields 



5X 



a x r A 4 + v x r A 



(3.28) 



where 5X denotes the deformation of the metric component, ax and vx are two constants and A is 
the dimension of the corresponding operator. In the dual theory, such a deformation corresponds to 
deforming the Lagrangian by £deform = £cft + a x@, where O denotes the corresponding operator 
and vx = (O) denotes the vacuum expectation value (VEV) of the field. 



From our solution, especially using the expressions in (3.26) and (3.27), it can be shown that 



the deformed metric functions have the following asymptotic behaviors 



(3.29) 



where /3f tg is some dimensionless numerical constant that will not be important for us, and r* is 
the UV cut-off. It is clear from the above expansion that m -2 corresponds to a coupling which 
has dimension —2 in the dual field theory. Thus the deformation to the CFT in this case will be 
of the form: £deform = £cft m 2 ^> where dim[C] = 6. The non-normalizable term in the above 
expansion of the metric function means we have inserted an irrelevant operator of dimension 6 
in the dual field theory. This fact conforms to the conventional wisdom that flavor back-reaction 
generally gives rise to irrelevant deformations of the background [27J. 

Let us be more precise about the identification of the various fields in the gravity background 
with the operators in the dual field theory^] Following the discussion in section 3.4 of [27], the 
scalar fields that are dual to the dimension 6 and the dimension 8 operators — which are denoted 
respectively by p and q — are given by 



1 



P 



(f-g) 



q 



2 / 5 

-lf + 4g +l \0 g h 



(3.30) 



Now using the most general solution that is presented in (|C.44|)-(|C.48|), we get 

2 



P 



e — r + 
m 2 



g=-logL + -| — + 



3 VL 4 



hr 4 H 



rl + • • . 



(3.31) 



The choice made in (C.49) is equivalent to setting the source for the dimension 8 operator to zero. 



Note that, as far as the gravitational background is considered, the deformations are very much 
similar to the ones discussed in e.g. [27]: i.e. because of the running dilaton, we have a deformation 
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We are thankful to Javier Tarrio for pointing this out to us. 
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from the conformal Klebanov-Witten solution to the non-conformal one. There is also a relative 
metric deformation between the S* 2 x S 2 base and the U(l) fiber direction of the original T 1,1 . As 
is demonstrated in e.g. [10], the departure from the conformal invariance is sourced by TrF 4 -term 
in the dual CFT, which is a dimension 8 operator. As for the identification of the dimension 6 
operator (discussed in the previous paragraph) is concerned, we will remain agnostic since the 
Lagrangian of the probe sector itself is not well- understood at present. However, this will not 
prevent us from exploring the phase diagram of this theory as the coupling m -2 is varied, 14 



3.2 A few comments on the regime of validity 



To comment on the regime of validity of the perturbative solution, let us begin by recalling the 
solution for the dilaton 



log 



(3.32) 



Quite clearly, the dilaton blows up exactly in the limit r — > oo. Thus the dilaton in our first order 
back-reacted background diverges only at infinity, which means there exists a Landau pole when 
the UV is located strictly at infinity^ This is similar to the situation studied in [32] where the 
Landau pole appears at r — > oo if we are to take the leading order solution in e. Thus, as far as 
the Landau pole is concerned, we can safely use our solutions for r* <C oo, where r* is the UV 
cut-off. 



Now to make sense of the perturbative solution given in (3.23)-(3.27) for the range of radial 
coordinate r# < r < r* we must impose 



log 



m 



r_H_ 



< 0(1) and 



li 



< 1 



e lm~ 2 r 2 | < 1 



(3.33) 
(3.34) 



along with the condition that e < 1. The first condition comes from demanding that the pertur- 
bative solution for the dilaton makes sense and the second two conditions come from demanding 
the same for the functions e* and e 9 . 

As discussed in [32], such back-reacted solutions can have other UV pathologies. Such patholo- 
gies show up in the so-called "holographic a-function" . It can also be checked explicitly that the 
holographic a-function in our case always remains finite and there is no additional UV cut-off in 

14 Perhaps a better expression would be "the family of theories" , since varying a coupling by hand means we are 
continually changing the theory. 

15 We will elaborate on this point in the next section. 
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the game. This is most likely because we found a solution only up to the leading order in the 
back-reaction. Thus we have a simple hierarchy of scales in our system 



r H < r* < r LP 



with r L p — > oo , 



(3.35) 



where the subscript "LP" stands for Landau pole. 

The validity of the SUGRA+DBI etc will give identical conditions as in [32]. Suppression of 
closed string loops imposes N c 3> 1 and the suppression of a' corrections imposes A ^> 1. For the 
smearing procedure we also need a large number of probes, which means Nf ^> 1. We will hold 
Nf/N c fixed and this is essentially the parameter that we control. Also, we need to impose that a' 
corrections to the supergravity background (which comes from leading order terms like a /3 i? 4 and 
hence scales as A _3//2 ) are sub-leading compared to the flavor back- reaction. This further imposes 
that \- 3/2 < e. 

3.3 A few comments on the dual field theory 

In this section we will offer some comments regarding the dual field theory. We will closely follow 
the discussion presented in [27]. Before taking the back-reaction into account, the field theory 
dual to the AdS 5 x T 1,1 is an = 1 superconformal quiver gauge theory with gauge group 
SU(N C ) x SU(N C ). This theory has two bi-fundamental matter fields, usually denoted by A4 and 
Bi, i — 1, 2 with a quartic superpotentialjS] 



where A is the 't Hooft coupling. This theory has a global SU (2) x SU (2) x U (1)r x Z2 symmetry. 

Adding the probe sector a la [20J introduces chiral flavors and realizes a global U (Nf) l x U{Nf)n 
flavor symmetry, where Nf is the number of D7 and anti-D7 branes. At zero temperature this 
flavor symmetry is spontaneously broken down to a U(Nf)d ilLg ^IU\, but at finite temperature this 
symmetry is restored (38]. Since this theory is not supersymmetric by construction, it is not 
currently clear how the Lagrangian of the probe sector should look like. For some thoughts 
regarding this issue, see [20]. Nonetheless, a lot of interesting physics can be extracted purely 
using the gauge-gravity duality. 

When considering the back-reaction, we are modifying the theory. The general procedure of 
taking into account back-reaction by the flavor sector introduces localized sources and therefore 
the global symmetry is expected to be of the form U(Nf)i x U(Nf)n. However, to avoid the 
complexity of the problem and to have good control we have taken advantage of the "smearing 
procedure" as outlined earlier and in appendix A. The smearing procedure distributes the branes 
and the anti-branes such that we recover the full global symmetry of the internal manifold that 



W = \Ti(A l B 3 A k B l )e ik e= l 



(3.36) 
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we had before introducing any flavors. The smearing is done homogeneously such that there is 
one D7 and one anti-D7 at each point of the two S 2, s. Thus the flavor symmetry group is broken 
down to two copies of 



U(N f ) L x U(N f ) R -> U(lf L f x U(l)% 



(3.37) 



Now, we are considering the back-reaction at finite temperature. The presence of a non-zero 
temperature breaks the conformal symmetry in the sense that it introduces a scale in the system. 
In the probe limit, the flavor sector breaks this conformal symmetry spontaneously even at zero 
temperature [2D] . Thus the inclusion of the back-reaction will explicitly break conformal invariance 
and will lead to a non-zero beta function. In order to make this precise, we need to know the 
relations between the gauge couplings, the theta angles and the supergravity fields. Such relations 
can be derived properly for the M = 2 orbifold theory which is dual to the background obtained 
by placing a large number of D3-branes at the singularity of C x C 2 /Z 2 . In this orbifold theory 
such relations are given by 



4tt 2 
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7re 
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~9l~ 
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9s 
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-irC 


1 

~ 27 


I c - 

Js 2 
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B2 — n (mod 2tt) 
(mod 27r) , 
(mod 2tt) , 



(3.38) 
(3.39) 
(3.40) 
(3.41) 



where gi^ and are the two gauge couplings and theta angles corresponding to the two gauge 
groups, B 2 is the NS-NS two-form, Ci is the RR two-form and Cq is the axion field. Here the 
integrals are performed over the two-sphere which shrink at the orbifold fixed point. We will 
assume that the above relations hold true for the conifold theory as well,^ 

In our back-reacted background, we do not have any if 3 = (IB2, C2 or Cq. Assuming that the 



two gauge couplings are equal (g 



YM 



#i7 2 = #1/2), we get 
4vr 2 



9ym 



71 



(3.42) 



Using the definition of the 't Hooft coupling A = gy M N c , the solution for the dilaton from equation 



(3.24), the definition of e from (3.22) and the above relation we can compute the beta function to 



16 Howcvcr, in 44 it has been suggested that the formulae relating the sum of the gauge couplings and the 
theta angles should be modified. Nonetheless, in the spirit of [57] we will use these formulae to extract qualitative 
features. 
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be given by 

r/3 A = r (d r X) = (J^j ^A 2 > . (3.43) 

Thus the beta function acquires a positive contribution of the order Nf/N c , which implies that 
there is a Landau pole in the UV and thus this theory lacks an UV completion. The location of 
the Landau pole is exactly where the dilaton diverges, i.e. strictly in the r — > oo limit. 



4 Thermodynamics 



In this section we will briefly discuss the thermodynamics of the finite temperature back-reacted 



background given in equation (3.14) and (3.23 )-(3.27). To make connection with earlier works 



we will keep our discussion closely analogous to what is discussed in e.g. [32]. The Hawking 
temperature of the black hole, which is identified with the temperature of the dual field theory, 
can be obtained by first compactifying the Euclideanized time-direction and then requiring the 
regularity of the near-horizon metric. The temperature is identified with the inverse period of the 
compact Euclidean time-direction. Thus we get 



T 



1 d 



(h,- 1/2 (r)b(r)) 



L 2 



7T 



47T dr 

The entropy of the ten-dimensional black hole is given by 

2ttA 8 27 



1 - 



"10 



32 



-7T 2 iV 2 T 3 



1 

1 + 



16 



£>(e 2 



(4.44) 



(4.45) 



where Ag denotes the eight-dimensional area of the even-horizon transverse to the time and the 
radial coordinate. In obtaining the above result we have also used the fact that 2k 2 = (27r)V 4 ^ 2 . 
Note that, the above result is not completely new in the sense that we could obtain the same 
by computing the entropy of the background (before taking any back-reaction) and the probe 
separately and adding them up. 

Proceeding along the same direction, we can evaluate the ADM energy of the black hole 
spacetime using the formula 



«10 ^RS J 



(4.46) 



where G tt is the time component of the background metric, G% denotes the metric on the constant 
time and constant radius hypersurface and finally Kt and Kq denote the extrinsic curvatures of 
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the eight-dimensional hypersurface within the nine-dimensional constant time hypersurface at non- 
zero and vanishing temperatures respectively. The formula for evaluating the extrinsic curvature 
is given by 

K:= -j=L= d^y/fetGgn^ , where n" = -^=6? , (4.47) 

where Gg denotes the metric on the constant time hypersurface. 

The subtraction of the zero temperature background should be viewed as a regularization 
scheme. Note that, in our case, the zero temperature back-reacted solution is likely to be more 
complicated, so we have not pursued it here. However, we are going to use the above formula for 
computing the ADM energy and we will simply send th — > while computing Kq. Operationally 
this will simply get rid of the (quartic) divergence term and thus will amount to effectively reg- 
ularizing the energy function. Moreover — even though the back-reaction at zero temperature 
potentially convolutes the analysis — far in the UV the probes have qualitatively similar behavior 
both at zero and at finite temperature. Thus, purely on physical grounds, it also makes sense to 
subtract the term corresponding to Kq which is evaluated by sending r# — >■ 0. With all these 
prerequisites, the final result is obtained to be 

S1 ^N 2 c T A (l + ]-e)+0(e 2 ) . (4.48) 



64 

We observe that the correction due to the presence of the flavor degrees of freedom contributes 
equally to the entropy and the internal energy. With this result, we can use the standard definitions 
and rules of thermodynamics to evaluate e.g. the specific heat, cy := 8t£ or the Helmholtz free 
energy, F = e — sT. 



5 Probe branes and phase transition 

In this section we will discuss the physics of an additional probe sector, a la mode [20], that can 
be introduced in the back-reacted background given in (3.23)-(3.27). Our task here is to analyze 
how the back-reaction changes the physics in this additional probe sector. 



5.1 D7 and anti-D7 branes: phase transition 

We will add Nj additional D7 and anti-D7 branes in this backreacted geometry with N'j- <C Nf and 
we will ignore the backreaction of these additional branes. We will study a D7-brane configuration 
which spans the space-time coordinates x M , the radial direction r, and the coordinates ?/>,#2,02- 
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The transverse directions are given by two-sphere coordinates Q\,<\>\. The DBI (Euclidean) action 
is given by 

•S'dbi = t 7 [ e^^g~ 8 d 8 x . (5.49) 



Here, g% is the determinant of the induced metric. Assuming Q\ = 6(r), <pi = 0(r) we get 

S DBl = M T J dre sir)+29{r)+ ^ r) ^jh(r)e 2 ^) (6' 2 + sin 2 9^) + 1 , (5.50) 

where Mt = i^N'^^V^i^Ti 2 ) / '(9T) and T is the background temperature. The factor of 2 appears 
because the brane and the anti-brane contribute equally to the thermodynamic free energy. The 
equations of motion obtained from the above action admit two classes of solutions: Firstly, the 
parallel-shaped solution 

cj)(r) = const, 0(r) = tt/2 , (5.51) 
and, secondly, the U-shaped solution 

6(r) = | , (5.52) 

6c 

(h'(r) = — (5 53) 

V ' v /6( r )2 e 2/(r)+8 5 (r)+2*(r) _ Q c 2^ r y2g{r) ' V ' ' 

where c is some constant which can be determined by using the condition that 4>'(r) diverges at 
some radial point r = r 

The asymptotic behavior of the profile function is given by 

0(r) = ^-^ + ... . (5.55) 

* 

Thus the asymptotic angle separation is the non-normalizable mode (corresponding to a source) 
and the constant c is the normalizable mode (corresponding to a VEV). The qualitative picture of 
the possible embedding functions is again similar to the one demonstrated in fig. [TJ The parallel- 
shaped represents a chiral symmetry restored phase and the U-shaped ones represent the symmetry 
broken phase. 

Before proceeding any further, we can explore how the coupling A0* depends on the various 
parameters (e.g. e or m) in this back-reacted background. A representative plot is shown in fig. [2] 
A couple of comments are in order: the constant c grows monotonically with r and therefore 
the dependence of A0* vs c will be qualitatively similar]^ We have chosen to plot r since the 
corresponding values for c are very large numbers. 



17 Note that in the thermodynamic ensemble, A<^* and c are canonically conjugate variables and therefore we 
should be able to conclude about the existence of a phase transition by looking at this plane. 
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Figure 2: The dependence of the asymptotic angle separation with the radial variable ro where 
the brane-anti-brane pair meets. The blue and red curve corresponds to m~ 2 = 0.06 and —0.06 
respectively. Here we have set r# = 1, r* = 10 and e = 0.01. 

In [38] it was demonstrated and subsequently argued that in the original Klebanov-Witten 
background at finite temperature, the coupling A</>* monotonically increases with increasing ro 
and approaches a constant non-zero value where the curve flattens. This in turn means that there 
can be no phase transition. This is in stark contrast with what we observe here: From fig. [2] 
we see that the A0* curve tends to flatten for r considerably larger than the horizon radius; 
however, then the curve starts decreasing as ro becomes closer to r* and eventually goes to zero. 
This signals that there will be a phase transition once the back-reaction is taken into account. We 
also observe that for m -2 = 0.06 the curve bends faster than for m~ 2 = —0.06; thus we expect 
that the critical A0* — separating the two phases — will increase as we go from m -2 = 0.06 to 
m~ 2 = —0.06. Note that this qualitative difference in behavior is sourced by the fact that we need 
to use a finite cut-off r* and therefore it is the presence of the running dilaton in our back-reacted 
background which changes the physics drastically. 

We can also learn about the existence of possible phases from this curve: this is schematically 
shown in fig. [3} Clearly, there exists a maximum A0* for a given m beyond which only chiral 
symmetry restored phases are available. On the branch where (<9A0*)/(<9ro) > 0, the system is 
thermodynamically unstable. This is because increasing ro by a small amount will increase the 
angle separation and will either push the brane-anti-brane pair to infinity or pull them all the way 



19 



A0„ 



2.0 


*SR only 


1.5 




^SB possible 


1.0 


^SB possible 




0.5 








2 4 6 8 


10 12 



^0 



Figure 3: The dependence of the asymptotic angle separation with the radial variable ro where 
the brane-anti-brane pair meets. Here we have set r# = 1, r* = 10 and e = 0.01 and m~ 2 = 0.06. 
The horizontal dashed black line represents the maximum value A</>* beyond which there are no 
chiral symmetry broken phase. 



to the horizon. However, note that unlike the case encountered in [3B], chiral symmetry breaking 
is possible for the entire range < A0* < A0* max , since the A0* curve bends down all the way 
to zero. 



Now, using the equations of motion we can calculate the on-shell action,^ For the first set of 
solution, we have 

poo 

dre f{r)+2g(r)+*(r) _ ( 5 5 g) 



S\\ = M T 

For the U-shaped solutions, we get 

Su=H T 



f rfre /(r)+29(r)+<I>(r) ^6(r)e 2 fM0'(r)2 + 1 



(5.57) 



with 4>'{r) given by equation (5.53). To find out the energetically favored embedding, we need to 
look at the difference in free energy of these two classes of embeddings: 

AS = {S v - S\\)/M T . (5.58) 

To find out the phase diagram, we will evaluate this quantity numerically. 



18 Rccall that in the Euclidean signature, the on-shell action corresponds to the thermodynamic free energy of 
the corresponding phase. 
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To represent the phase diagram, let us first define the following dimensionless coupling 

™ 2 = ^ i 1 + k) ■ (5 ' 59) 

Thus the coupling m 2 measures the relative strength between the background temperature and 
the dimensionful coupling m. We will represent the corresponding phase diagram in the A0* vs 
m 2 plane. The resulting phase diagram is shown in fig. [4] and fig. [Hj 



A0, 




Figure 4: We have shown the phase diagram in the Ac/)* vs m 2 plane. The blue solid line separates 
the xSB = chiral symmetry broken broken phase (below the line) from the xSR phase (above 
the line). Above the red solid line only xSR = chiral symmetry restored phase is available. We 
have shown these in more details in fig. [5j We have used e = 0.01, rjj = 1, L = 1 and r* = 10. 

From fig. [4] we observe that increasing m 2 , and therefore increasing temperature with respect 
to m 2 , favors the chiral symmetry restored phase. This is expected since the introduction of the 
coupling m~ 2 induced a spontaneous breaking of chiral symmetry which would otherwise remain 
restored at finite temperature. Thus the irrelevant operator that we have introduced via taking the 
leading order back-reaction into account has induced this spontaneous symmetry breaking. The 
detailed phase structure is magnified in fig. |5j From this we also notice that — comparing with the 
corresponding phase diagram that we analyzed in [38] — there is no "runaway or metastable x^R" 
phase available when the back-reaction is non-vanishing. Thus the back-reaction "stabilizes" a 
part of the phase diagram and makes it physically accessible. 
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A0. 




Figure 5: We have shown the details of a part of the phase diagram presented in fig. |4| Above 
the red dashed line we only have the xSR phase; between the red dashed we have both xSB and 
%SR as available phases. The blue solid line separates these phases. We have used e = 0.01, 
ru = 1, L = 1 and r* = 10. 

If we are to dig deeper into the origin of the non-trivial phase structure that we obtained here 
we may be tempted to conclude the following: taking the back-reaction into account has explicitly 
broken the conformal invariance and therefore we can construct a meaningful quantity, e.g. the 
ratio T 2 /m 2 , which can favor one possible phase over the other depending on its strength. This is 
indeed true; however, it is not the complete story. 

Recall that the back-reaction considered in e.g. [32] - which is different from what we have 
considered here — also introduces an irrelevant deformation to the original CFT and introduces 
a non-trivial scale. Forgetting about the possible interpretation for the moment, we can imagine 
introducing an additional set of ATJ. probe D7 and anti-D7 branes in the finite temperature back- 
reacted background obtained in [32J . However, it is straightforward to check that in that case, the 
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so-called phase always remains favored and a phase transition never happens. This indicates 
the phase transition we observe here has more to do with the details of the theory, i.e. precisely 
what irrelevant operator we introduce, than just an effect induced by the breaking of the conformal 
invariance. 

5.2 Thermodynamics of probe D7-branes 

After discussing the physics of the phase transition, we will now discuss the associated thermody- 
namics in this probe sector. 

5.2.1 Parallel-shaped embedding 

To begin with, we will analyze the parallel-shaped (also sometimes called the black-hole) embed- 
ding. In this case, it is possible to obtain analytical results and hence we will discuss this case in 
detail and then move on to the discussion for the U-shaped profiles. The profile is described by 



const, 



0(r) = vr/2 . 



(5.60) 



The Helmholtz free energy F is given by TSu, where Su is the on-shell Euclidean action for the 
parallel embeddings: 



(5-61) 



The overall constant Mr has been defined earlier, e.g. below equation (5.50). The on-shell action 
diverges as the regulator, denoted by A above, is taken to a large value 
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(5.62) 



Setting the UV-cutoff A = r*, we get 
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(5.63) 



Apparently, there are additional and r 2 divergences along with the usual quartic divergence. 
However, we have to be careful in treating "diverging" terms which are multiplied by our expansion 
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parameter e. In light of the discussion in (3.33)-(3.34), it is clear that we only have a quartic 
divergence rj We can remove the quartic divergence by adding the following counter-term: 



5 ct = (2N' f r 7 ) 



-1 



1 1 2 -2 

32^3^ 



d 7 x\J det [7] 



(5.64) 



where 7 denotes the induced metric on the probe at a constant radial slice. So, finally we have 
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The background temperature is given by 



rp r H ( 1 

Therefore, using the definition of the 't Hooft coupling A = 4ng s N c , we get 

1 



(5.65) 



(5.66) 
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(5.67) 



Above, we have defined So = \ [ ). Now we define A e ff(T) 



Aeff(T) = A 



1 + i el0g (^) + k V * e 



+ 0(6 5 



(5.68) 



This defines for us an effective 't Hooft coupling at finite temperature. Using this effective 't Hooft 
coupling we get, 

5|, + 5 ct = -s {N>V R3 T 3 \l G ) (l + ^ + 0(e 2 ) . 
The free-energy density can be obtained from the on-shell action 



(5.69) 
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-s N' f T*\ cS (T) 2 (l + ^e)+O(e 2 ^ 



Consequently the entropy density, energy density and pressure are given by 



H = 4s N>T 3 \ cS (T) 2 (l + lie ) + 0(e^ 



e„ =/), + T S]l = 3s N' f T*\ cS (T) 2 (l + ^e) + 0(e 2 ) 

Pi, = - /I, = , iv;T 4 A eff (r) 2 (l + A e ) + 0(e 2 ) . 



(5.70) 

(5.71) 
(5.72) 
(5.73) 



5 We will deal with the logarithmic term momentarily. 
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5.2.2 U-shaped embedding 

Next, we will discuss the thermodynamics for the U-shaped embedding described by 

0(r) = \ , (5.74) 

fir 

Aj( r ) — C5 75) 

V ' v /5( r )2 e 2/(r)+8 9 (r)+23>(r) _ Q c 2^ r y2g{r) ' K ' 1 



where the constant c is determined via equation (|5.54|). The on-shell action is given by 

Su=N T 



dre f(r)+2 9 (r)+*(r) Jl&( r ) e 2ff(r)0/( r )2 + 1 _ (5 ?6 ) 



Again all divergences can be removed by adding the counterterm: 



S ct = (2N' f )r 7 j (-1 + l e - l^m- 2 ) J d 7 x^d^r\ . (5.77) 



Now all the thermodynamic quantities can be computed: 



h =f\\ + (^Ti X2N 'f) 5 ^ ' ( 5 ' 78 ) 
*u =*|| + f-lnA 2 ^ Ss(T) , (5.79) 



y 21 

eu =e„ + (^|nA 2 iV;) 5e(T) . (5.80) 

Here 8f(T), Ss(T) and Se(T) needs to be evaluated numerically and are shown in figs. [6] From 
fig. [6] we observe that the free energy difference crosses zero at a certain value of the temperature 
measured in units of r , which encodes the existence of the phase transition. We also observe that 
the entropy is always greater in the chiral symmetry restored phase. This is also expected: in 
the chiral symmetry restored phase more degrees of freedom are available since they cannot be 
bound to form a chiral condensate. This fact also gives rise to an increased internal energy for the 
symmetry restored phase. 

A first order phase transition is accompanied by a jump in the entropy density at the critical 
temperature that results in a non-zero latent heat. This latent heat is simply given by 

-^latent = T c 5s(T c ) . (5.81) 

On dimensional grounds, the overall scaling of the latent heat should behave as T c 3 . However, here 
we have other scales in the system and thus the most general behavior the latent heat can follow 
will be of the following form 

(T 2 T \ 

L latent =T c 3 F( ^,^J , (5.82) 
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{100 Sf, 5s, 10 Se} 




Figure 6: Temperature dependence of 1 00 5f (T) (black), 5s(T) (blue) and 10<5e(T)(red) for e = 0.01 
and m -2 = 0.01 and L = 1. Along the horizontal axis the temperature is measured in units of 
No- 
where F is some unknown function and A is the U V cut-off of our theory. In fig. [7] we have shown 
a representative behavior of the latent heat with the critical temperature. 



6 Some comments about the stability of the additional 
probes 

Since we are dealing with non-supersymmetric probes to begin with, we would like to understand 
whether after taking into account the effect of the back-reaction, this system remains stable or not. 
Ideally speaking, one needs to study the fluctuations of the full background to make any concrete 
claims about the stability. We will, however, do the following: we assume that any possible 
instability in this system will be sourced by the non-supersymmetric probes and by analyzing the 
stability of the additional set of probe branes in the back-reacted background, we should be able 
to see the effect of this instability, if there is any. Here we will perform suggestive calculations to 
shed light on the issue of stability of the additional probe branes. 
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Figure 7: Temperature dependence of the latent heat during the phase transition. Here we have 
used e = 0.01, m~ 2 = 0.01, L = 1 and we have also set ro = 1. 



6.1 Computing the force on the branes 

One way to ensure stability of probe branes is to simply embed them in a supersymmetric way, 
which is not the case we currently have. One easy way to see e.g. a supersymmetric "extra" set of 
probe Dp-branes is stable in the background of a large number of branes of the same dimension is 
to compute the probe action assuming that the probe brane is "moving slowly" in the transverse 
directions. In this case, a cancellation occurs between the Dirac-Born-Infeld piece and the Chern- 
Simons piece and the resulting Lagrangian takes the form of that of a free particle at the quadratic 
order in velocity. 

We could ask the same question in our non-supersymmetric set-up as well. What we will try 
to explore here is what happens if we give this additional set of probes some velocity along its 
transverse directions. In our back-reacted background, we do have a "large" number of D7-branes, 
but we do not have any non-trivial RR-potential sourced by them since there is an equal number 
of D7 and anti-D7 branes uniformly smeared. Thus the additional probe D7-brane will not have 
any contribution from the Chern-Simons term and it's action is entirely given by the DBI piece. 
Thus, even in the small velocity approximation, the effective Lagrangian will not be that of a free 
particle - but will be of a particle moving in a non-trivial potential. Our goal here is to determine 
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this potential. 

To make this more precise let us consider, in the back-reacted background, an additional probe 
D7-brane whose embedding is given by 

9 = 9(r) + vgt , = 0(r) + V . (6.83) 

Thus we are considering a very general embedding where both 9 and <fi can be a function of the 
radial coordinate. We will try to extract the Lagrangian in the limit when the velocities are small. 
Note that, in the background there is an U(l) isometry along the direction 0. Thus without any 
loss of generality we can set = 0. 

Now, the Lagrangian that results from the above ansatz is complicated looking. However, 
one can check the following: in the limit of small velocities, the leading order term (which is 
independent of the velocities) determine the classical profile of the probe and not surprisingly a 
simple solution for the classical probe is given by 

9 {0) (r) = n/2, (o) (r) = const . (6.84) 

Here the subscript stands for the classical solution. To have a systematic expansion, in the small 
velocity limit, the profiles will deform and we denote this by 

9 = 7r/2 + a (6*(i)(r) + vq£) , = const + a0(i)(r) , (6.85) 

where a is a small parameter. Now it can be explicitly checked that the DBI Lagrangian takes 
the following form 

£ = £ (o) (r) + £ (1) (^,^ (1) ,0 / (1) ,r)a 2 + O(« 4 ) . (6.86) 

Thus the equations of motion obtained from the Lagrangian at the order a 2 gives the solution 

0(1) = 0, 0(d = O. (6.87) 

Furthermore, there is no non-trivial potential (which depends on the direction 9) and thus the 
velocity vq remains a constant of motion. Thus, a small velocity along either the <fi or the 9 
direction does not get accelerated suggesting that the embeddings are possibly stable. 

It turns out that one can generalize this argument to a full non-linear level. With the ansatz 



in (6.83), the DBI Lagrangian is given by 

e %0'2 + i (66 - e 29 hvl) (6 + be 29 sin 2 9<P' 2 ) 
£(9,9',9,(j)') , 9 = vg. (6.88) 



r 1 p /+2g+<E 

~ 108 
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The equation of motion will now imply 

dC 



c = const . (6.89) 



Note that unlike the small vg case, this constant even for the probes that go all the way to 
the horizon q 

The other equation of motion will be given by 

d fdC\ d fdC\ dC 



Using the explicit form of the Lagrangian in (6.88), it can be shown that we still have the following 
solution 

9{r) = tt/2 , ^v e = 0. (6.91) 

Hence the static probes, when given an initial velocity along the angular directions, will not 
accelerate. This also suggests that the probe sector is indeed stable. 

6.2 Asymptotic behavior of the profile functions 

So far we have focussed on the solution 9 = n/2 and <ft = const. But we can also look for non- 
constant solutions and from their asymptotic fall-off behavior we can read off the dimension of the 
corresponding operator. Typically in AdS/CFT, an instability (due to non-unitarity) of the dual 
gauge theory is reflected in the fact that the dimension of this operator becomes complex-valued. 

Although the back- reacted solution does not asymptotically approach AdSs x T 1 * 1 , the AdSs 
part is still present. Thus we can use the above intuition to explore whether a non-constant 
solution for either 9 or has a problematic UV behavior. 

The equation of motion for <fi is given by 

e /+4g+<I> 6sin 2 qjj 



[l + {l/6)be 2 9 (9' 2 + sin 2 9<f)' 2 )] 1/2 
The equation of motion for 9 is given by 



C<s> . (6.92) 



d 
dr 



2^, 2 MV2 



[1 + (1/Q)be 2 a (9' 2 + sin 2 9(j)' 2 )] 1/2 



(6.93) 



[l + (l/6)6e 2 f (9' 2 + sin 2 0<jP)] 

20 Note that the solution corresponding to c = will run into trouble in the deep IR. It represents a brane that 
is moving with some velocity along the (9-direction at the asymptotic boundary. This velocity will be red-shifted 
as we go closer to the horizon and eventually exceed the speed of light. To prevent this, the brane will bend and 
develop a non-trivial profile e.g. along the <p direction. Thus <j>(r) will not remain a constant all the way. Such an 
effect has been discussed in e.g. |46j . However, in the small vg expansion this bending effect can be neglected and 
we can approximate the brane having a straight shape. 
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Now let us first check the case: 9 = it/ 2 and <fi = const + <f>(r). Assuming that (/)' — > fast 
enough as r — > oo such that we can ignore the (//-dependent term under the square root in the 
denominator, we get 

8( ' ; (6.94) 



r 5 (8 — e + 2elogr) 

Using a small e expansion, the asymptotic behavior can now be obtained to be 

4e log r — 16 — e 

<p = const + cs — — h • • • . (6.95) 

64r 4 

Thus the asymptotic behavior picks up a logarithmic term, but it does not have any monomial in 
r whose exponent is complex-valued. 

The other simplifying case we can explore is to consider 9 = n/2 + 9{r) and <fi = const. In this 
case also one can show 

4e log r — 16 — e , 
9 = n/2 + c e + ... . (6.96) 

Actually we can simultaneously consider 9 = n/2 + 9(r) and <fi = const + 0(r) and using the 
equations of motion it is straightforward to argue that their UV behavior are exactly the same as 
we have written here. The above argument is fairly suggestive that the probe should be stable. 

6.3 Studying the fluctuations 



Let us now study the fluctuation equations obtained in (|D.24|) and (|D.25|) 

d, 



e^V^gf^dbSO^ + e $ v /Z ^«9^50 1 = , (6.97) 
e^V^g^Gl^d^} - e*V=E*di,60i = . (6.98) 



Here G° denotes the background metric evaluated at the classical profile of the probe, go denotes 
the induced metric of the classical profile, E° denotes the on-shell determinant of the DBI-part of 
the action and the indices a, b run over the worldvolume coordinates of the probe brane. 

It is straightforward to check that for the given background, both the fluctuations obey the 
same equations of motion. We will henceforth denote this generic fluctuation mode by 5X. It is 
difficult to solve the above two equations in complete generality. It is also not possible to have 
separation of variables in the most general case. Thus we need to make some simplifying assump- 
tions to make any progress. 

Case 1. We begin by assuming 

d^SX = , (6.99) 
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where SX represents any of the fluctuation modes. Note that this simplification should not com- 
promise the generality of our discussion. Even after introducing the additional probe sector, the 
symmetry along the ^-direction is not affected and thus we do not expect any instability to arise 
along that direction. 

Furthermore, we assume that the fluctuations are independent of the Minkowski directions. 
With this simplification we can separate the variables 2 ^ and the resulting equations are obtained 
to be 



d r [e f+i9+ ^bd r Y{r)] - 6£(£ + l)e /+29+ *F(r) = 

2 . 
S 2< 



V% 2 Z + £(£+l)Z 



. 



(6.100) 
(6.101) 



Here we have used 5X = Y(r)Z(02,(f)2)- 

We can recast the radial equation of motion in the following form 




Or 




d T K(r 



where 



with 



And also, 



Y(r) 



6£(£+l)K{r) - V{r)K{r) 



a(r)n(r) 



a 
a 



a 2 , 1 ax 

— + -—dr 
a\ 2 <23 



ai(r) 
A(r) 

V{r) 



B(r) 



a-2 
a i 



a 2 (r) = d r (e /+49+$ 6) 

1/2 



1 

52 



A + 



03 

a 1 

B' 
~B 



-d r A 



~B 



A 
2 



A + 



~B 



(6.102) 

(6.103) 
(6.104) 

(6.105) 
(6.106) 

(6.107) 



Here ' = d r everywhere. 

The advantage of expressing the fluctuation equation in the form given in (6.102) is that it 
takes the form of a Schrodinger equation with a given potential denoted by V(r). By looking at the 
behavior of this effective potential one can draw conclusions about the nature of the fluctuations. 
Now we can look at how the effective potential behaves as we introduce the parameter e and the 
constant m, 22 This is shown in fig. km From fig. M it is clear that for m -2 < the effective potential 



21 If we assume a plane- wave type ansatz for oscillations along the Minkowski directions, then we can also separate 
the variables. 



2:2 



Note that in this case these two parameters appear in a single combination (em r H ) 
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V(r) V(r) 




(a) (b) 

Figure 8: We fix e = 0.01. Left panel: The behavior of the potential for various positive values 
of m~ 2 . The black, blue and red curves correspond to m -2 = 0,0.03,0.06 respectively. Right 
panel: The behavior of the potential for various negative values of m -2 . The black, blue and red 
curves correspond to m -2 = 0, —0.03, —0.06 respectively. 



tends to develop a runaway behavior for large r. For m -2 > 0, however, the effective potential does 
not develop any possibly pathological behavior as compared to the case of vanishing back-reaction. 
Actually it can be shown analytically that for large value of r, the effective potential behaves as 



V -> (4 + + 12rrT 2 er 2 + 0(r- 2 ) , (6.108) 

where r* is the UV cut-off. Now, recall that the condition that our perturbative solution is valid 
sets a constraint that e|m~ 2 r 2 | is parametrically much smaller compared to an order one number. 
This means that the effective potential should not become negative even for rrr 2 < 0. It is an in- 
teresting question how this effective potential might look like when e is not a small parameter. To 
answer this question, one needs to find the full solution of the equations of motion and this most 
likely calls for a numerical endeavor. Currently, we do not have enough information to comment 
on this situation. On the other hand, m~ 2 > clearly should not have any instability. 

Case 2: Let us now assume that the fluctuations are independent of the angular coordinates. 
Since there is an SO(3) symmetry, without any loss of generality, we can assume the fluctuations 
to depend on only one of the spatial directions. 
In this case, the equation of motion is given by 

d r (e f+i9+,s> bd r 5X) - ^e f+i9+ *d 2 t 5X + he s+A9+ *d 2 x 5X = . (6.109) 
Assuming a plane-wave type ansatz for the oscillations along the Minkowski directions the equation 



32 



of motion can be obtained to be 



d r (e f+A9+ Hd r Y(r-)) + u 2 ye f+4g+ *Y(r) - k 2 he f+49+ *Y(r) = , (6.110) 

b 



where 5X = Y(r)e~ iu)t+ikx . (6.111) 
6.3.1 Caseofw = 

We can again recast the radial equation in an effective Schrodinger equation with an effective 
potential. The process is very similar to the one outlined before, except the definition of the 
function a 3 will be different in this case, i.e. 

a 3 (r) = /ie /+49+ * . (6.112) 

The equation of motion now takes the form 

ai -d r ( x f^d r K(r)) - k 2 K(r) - V(r)K(r) = , (6.113) 



a 3 
where 

Y(r) = a{r)K(r) . (6.114) 
In this case the effective potential is given by 

V(r\ 6 ^ + 2 + rl ~ rfill-rt 

K > 32r 2 L 4 (r 4 - r A h ) 4r 2 L 4 (r 4 — rf t ) ' 1 1 

It is straightforward to check that this potential does not develop any qualitatively new behavior 
as e is dialed up. Of course, this is completely insensitive to the choice of the constant rrT 2 . 

6.3.2 Case of k = 

In this case, the function a 3 (r) is given by 

a 3 (r) = jV+ 4 s+* . (6.116) 
Here the equation of motion takes the form 




d r J— d r K(r) + u 2 K(r) - V(r)n(r) = , (6.117) 

«3 W a 3 J 

with the potential given by 

v( v _ 3 (r 4 - r 4 ) (3r 4 + 5r 4 ) e (-10r 4 r 4 + 9rj + r 8 ) 
{ ' 4r 6 L 4 32r 6 L 4 ' 1 ' 

Again, it can again be checked explicitly that this potential does not develop any new qualitative 
feature as we increase e. Thus there is no possibility of an instability within this analysis either. 
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7 Conclusions and future directions 



In this article we have studied the effects of back-reaction of flavor D7-/anti-D7-branes in the so- 
called Kuperstein-Sonnenschein model at finite temperature. We have treated the back-reaction 
perturbatively in which case it turns out to be possible to obtain an analytical solution for the 
background geometry. This can be regarded to be the first step towards a more involved compu- 
tation that will address the back-reaction non-perturbatively. Such an endeavor most likely calls 
for a numerical exploration where one needs to find series solutions for various unknown functions 
near the event-horizon and the asymptotic boundary and eventually construct a numerical solu- 
tion that successfully interpolates between these two limits. Such series solutions, especially near 
the asymptotic boundary, are rather non-trivial and take a general form discussed in [17] . It will 
be an interesting exercise to explore the numerical solution to learn more about the physics when 
Nf/N c is of arbitrary strength. This will shed light on the effect of the back-reaction for arbitrary 
values of Nf/N c on the phase diagram and the question of stability beyond the probe limit. 

Let us offer some comments regarding the stability issue. Before taking the back-reaction the 
naive picture for the construction is rather simple: we have localized probe branes and anti-branes 
located at two points in the 0-direction. They are equivalent to having Nf and — Nf charges at 
two points on a circle. A simple electrostatic intuition tells us these charges should attract each 
other and any little perturbation is likely to drive them towards each other and destabilize the 
system. However, our analysis does not confirm this intuition. Note that, contrary to what has 
been observed in [31] , we do not see any non-trivial cancellation between the DBI and the Chern- 
Simons terms in the quadratic fluctuation for the additional probe sector. Thus the stability that 
we observe is likely rooted in the fact that we have treated the back-reaction perturbatively. The 
mass spectrum of the system at finite temperature must be set by the scale of the background 
temperature. A deformation, brought about by a small perturbation, is unlikely to drive this mass 
spectrum to become tachyonic. This insinuates a limitation of our approach. 

On the other hand, the lack of a non-trivial cancellation as observed in [31] can be attributed 
to the use of the smearing technique which does not generate any non-trivial axion field to which 
the additional probes can couple to. Going beyond the smearing technique, one feasible scenario 
is that such a cancellation eventually guarantees stability for the system. Thus it is tempting to 
conjecture that in this system, in which supersymmetry has absolutely no role to play, various 
forces do perfectly balance each other to guarantee the desired stability. It will be interesting to 
understand whether such statements can be made more rigorous. 

As much as the above-mentioned directions of future work are promising to be very fruitful, 
there is more interesting physics to explore within the current setup. The additional probe sector 
introduced in the back-reacted background can teach us interesting lessons about this model. 
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In [38], we have explored the phase structure of the Kuperstein-Sonnenschein model at finite 
temperature and in the presence of an external electro-magnetic field. We can analyze similar 
physical questions, now including the effect of the back-reaction perturbatively. This will elucidate 
how the phase diagrams depend on the various parameters beyond probe limit computations. 
Another direction to pursue would be to include the effect of a chemical potential in this model. 
In the presence of such external parameters, we expect to recover a similar physical picture as 
reported in e.g. [3B1 HHl EH ED E21 E3]- Since the devoted reader has already experienced a fair 
share of agony and ecstasy, we will report on these in our forthcoming works [51] . 

Moreover, one could consider studying the back-reaction including the various gauge field 
excitations on the probes. Such gauge field excitations again correspond to a chemical potential 
or a constant electromagnetic field. Work along these directions has been performed in other 
models in e.g. [231 ESI E7J- It will be interesting to investigate and identify the robust universal 
features within our model as well. 

Note that the physics of chiral symmetry breaking can — perhaps more interestingly — also be 
realized by placing probe brane-anti-branes in the Klebanov-Strassler background [21] . However, 
the finite temperature version of the background is known only numerically in [UJ or approxi- 
mately for large temperatures [55]. The problem of considering back-reaction in this case will be 
technically more involved, but perhaps more rewarding from a physics point of view. Interestingly, 
in that model, the antipodal probe solution, simply described by <fi = const, exists even at zero 
temperature. Thus analyzing the back-reaction for such profiles may be a tractable and interesting 
problem to consider. 

Our work crucially depends on the smearing technique and it really captures the physics in an 
approximation — akin to the "s-wave" approximation — of the actual situation. Indeed we really 
have localized sources and not smeared sources. Although we realize interesting physics within this 
approximation, it is not evident how robust these features actually are beyond this framework, 
which is a technically more challenging problem because one has to consider localized sources 
and thus one will have to work with coupled PDEs instead of the coupled ordinary differential 
equations (ODE) that we have encountered herein. Nonetheless, it will be an interesting problem 
to consider in the futureF^I 

Let us conclude by saying that although the problem of understanding the strongly coupled 
physics of the quark-gluon plasma remains a daunting task, we are slowly making progress in 
uncovering some key universal features of the strong coupling regime using these model computa- 
tions. 

23 For some work from a slightly different perspective, where such localized sources have been considered, see 
e.g. [561 E3 [581159]. 
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Appendix A. The smearing form 

The idea of the smearing is to recover the full symmetry of the background we had before placing 
any probes. The Z 2 -symmetry can be recovered simply by picking another pair of D7/anti-D7 
branes with transverse plane described by the {# 2 , 02 } directions for each pair of D7/anti-D7 brane 
with transverse plane described by {#i,0i}. The induced line element on the radial direction of 
the probe brane is simply 

dr 2 , (A.l) 

which is invariant (up to reparametrizations) under the shifts: #i i2 —> #1,2 + #0 and 0i i2 —> 0i,2 + 0o 
where #0 an d 0o are arbitrary constant numbers. Thus we can initially pick # 1)2 = 7r/2 and 
0i,2 = 0( r ) to focus on the "great circle". 

We can smear the branes on this great circle by using the shift symmetry in 0i )2 ; then we can 
rotate this great circle using the shift symmetry in 0i i2 . This process should then cover the entire 
S 2 and thus recover the full symmetry of the background. The density of the branes is then given 
by 

N f 

AWoo = -^ SiU0 OO , ( A -2) 

where 6^ and 0^ are the asymptotic values of #i )2 and 4>i,2 respectively. 

Now to determine the smearing form we follow [35]. The most general 2-form invariant under 
SU(2) x SU(2) x £7(1) x Z 2 is given by 

Q = ai(r) (e 1 A e 2 + e 1 A e 2 ) + a 2 (r) (e 1 A e 2 + e 1 A e 2 ) 

+ a 3 (r) (e^g 1 ) A (e 3 + e 3 ) + a 4 (r) (e 2 + e 2 ) A (e 3 + e 3 ) + tu , (A.3) 



^ + (^ 2 ) 2 + sin 2 ^ 2 (0' lj2 ) 2 
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where a,, i = 1, . . . , 4 are four hitherto undetermined functions. Here u is another 2-form invariant 
under the same aforementioned symmetry and is independent from the other terms. Now u can 
be fixed in the following by demanding that Q is closed, i.e. dQ = 0. After some algebra it is 
straightforward to show that Q takes the following form 

Nth) N'f(r) 

Q = /v ' (sin 6 1 d0 1 A dxj>x + sin 9 2 d0 2 A d(j) 2 ) — dr A (dip + cos M^i + cos 9 2 d(j> 2 ) , (A.4) 

Air Atx 

which matches exactly the result of j4S], including the overall normalization constant. 

Now, in the notations of [45] , the embedding functions of the probe branes are given by 

fi = 0i,2 - 0oo , / 2 = 0i,2 - 0(r) , (A.5) 

where 0oo is the asymptotic value of 01,2 (which is fixed to be 7r/2 for a single set of branes). Now 
Q can also be evaluated using the formula [45] 

n = I Po^J (0, - 0oo) 5 {<t>i - 0(r)) d0oo#oo (^ A dfi) . (A.6) 
i=i ^ 

From the above formula it is clear that N'Jr) = (since there is no 7"0 -component for example). 
Thus Nf(r) is simply constant. At non-zero temperature the energetically favored embeddings 
are the ones which go all the way into the horizon [2H] arid therefore Nf(r) = Nf for all values 
of r. On the other hand, at zero temperature we know that the probe banes meet at some radial 
position given by ro (we fix this parameter in our problem). So for a given tq, Q = for r < r$. 
Thus we conclude 

N f (r) = N f Vr > r , 

= Vr < r . (A.7) 



Appendix B. The equations of motion 

To set up the conventions, we will start with the action for type IIB supergravity together with the 
D7-brane DBI action written in the Einstein frame. In our framework, we have an equal number 
of D7 and anti-D7 branes uniformly smeared and therefore we do not source any C& potential (the 
Hodge dual of the axion field). The action reads 



S — SsiJGRA + ^flavour 
1 



2 K io 



d W X y/-G w 



Nt 



2r 7 J2 d^e 



-G 



(i) 



R- l -d M ^d M ^-\\F^ 



-G 



(2) 



(b.i; 
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Here kiq is related to the ten dimensional Newton's constant Gn via 



2k 2 10 = !6irG N = (2tt)V 4 ^ . (B.2) 

In the supergravity part of the action, R is the Ricci scalar, G\q is the background metric in Einstein 
frame, $ is the dilaton and the five form field strength F 5 satisfies the self-duality condition: 
*F 5 = F 5 . 

(1 2) 

In the DBI part of the action G 8 is the induced metric on the seven brane and we have 
defined g s = e**, where is the asymptotic value of the dilaton field ,_ With this definition we 
have 



T7 



-(2tt 



rV 4 



(B.3) 



such that 



1 



2nf 



9s 



The factor of 2 in front of the DBI piece comes from the fact that we are adding Nf D7 and anti-D7 
branes and they contribute equally to the action. The relative factors of g s can be understood as 
follows: the supergravity part of the action comes from the closed string sector and thus contains 
an overall factor of g~ 2 whereas the DBI piece emerges from the open string sector and therefore 
contains an overall factor of gj l . Note that here we will consider the back-reaction of the probes 
which are described by &i — tt/2 and <pi = const with i — 1 or 2. 

Applying the smearing procedure in the transverse directions leads to 



N f 

E / ^ 



-G 



(1,2) 



— >■ 



Nf 
Air 



d^xe* sin 6 1 a/ -G 
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(B.5) 



where G 8 ' represents the induced metric on the probe whose profile function is given by #2,1 = vr/2 
and 02,i = const. The resulting ten dimensional action is given by 



S 



1 



?K 2 



d w x J-G 



10 



R 



d M $d M $- 



2r 7 N f 
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2r 7 N f 
Air 



d 10 xe* sin^ 
d w xe^ sin 6 2 



-G 



(i) 



—G 



(2) 



The equations of motion read 



E 



MN 



D M 8 M <5> 



T 



(i) 



MN 



+ T (2) 
^ 1 MN 



,71(3) 
MN ' 



2k1 (2t 7 ) N, ^ 



sin^V-G^ + sme 2 \/-G 



.(2) 



dFz = 



(B.6) 

(B.7) 
(B.8) 
(B.9) 



4 This implies that we need to impose the boundary condition $ — > when we obtain the solutions. 
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where 



Emn — Rmn — ^RGmn , (B.10) 



T (i) 

1 MN 



^(d M ®d N <S> - l -G M Ndp<bd p §) , (B.ll) 

T (2) _ p p PQRS ' t? w PQRST r m -, 9 \ 

1 MN - gg MPQRSF N ~ 8 X 5! \ P Q RST ^MN, 

t ( 3 )mn = 2^^gflavor = _\ (2N f )g* y J^G&G^S? 5* ■ (B.13) 

We have now denoted the string coupling constant, g s = g*. The superscript will be a helpful 
bookkeeping device that there will be a cut-off radial position, denoted by r = r*, where the 
boundary theory will be defined in the back-reacted background. When the back-reaction vanishes, 
this cut-off surface r* can be taken to infinity. 

In order to solve the equations of motion, we will work with the following ansatz (in Einstein 
frame), which has the full SU(2) x SU(2) x U(l) x Z 2 invariance 

( J„2 2g(r) 

ds\ = h- l l\r) \-b{r)dt 2 + dXida?] + /* 1/2 (r) ^- + C — V [dd] + sin 2 0, # 2 ) 

w) 6 rrv 



2 



F 5 = k[r)h[rf^ (e* A A A e* 3 A e r + A e Sl A e* 1 A e 9a A e* 2 ) , (B.I5) 
where the vielbeins are defined as 

e t = h -V4 b W dt ? e - 1 = h-^dx* , e r = /i 1 / 4 ^ 1 / 2 ^ , (B.I6) 

e V = h 1/4 e f (dijj + cos + cos M0 2 ) , (B. 17) 
3 

e dl - 2 = ^h 1/4 e 9 d9 1>2 , e^ 2 = ^h 1/4 e 9 sin 9 h2 d(f) h2 . (B.I8) 

V6 V6 

Here the functions h(r), b(r), f(r), g(r) and $(r) are unknown functions which we have to obtain 
by solving the equations of motion. The Bianchi identity dF 5 = is trivially satisfied and the 
quantization condition 

f F p+2 = N c r p for p = 3 (B.I9) 



I 

2 K io 



yields 

27 

A;(r)/i(r)V 9(r)+/(r) = 27ng* s N c l 4 s = 4L 4 := -j-\l A s , (B.20) 
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where the 't Hooft coupling A is defined as A = 4irg*N c . In deriving the above relations we have 
also used the fact that 



— (2n)- p a'- {p+1)/2 , (B.21) 

9*8 



where a 1 is the string tension. Thus the function k(r) is determined in terms of h(r),g(r) and 
f(r). Moreover, 



G^' 2) = -4 sin (^,i)e /(r)+2fl(r) , (B.22) 



1 



^/ Z ^ = wV / %)sin^ 1 sin^ 2 e / M +4 ^ . (B.23) 
10o 



The resulting equations of motion are 



25 



$"(r) + $'(r) f/'(r) + Ag'(r) + = ^Me^" 2 ^ , (B.24) 

4 (fe'(r)^(r) + b(r) (/'(r)^(r) + /(r) + Ag'(r) 2 ) + 2e 2/(r) - 49(r) - Qe' 2 ^) + 1 /J ^ S = 

(B.25) 

b'(r)h'(r) 16^8 -2(/(r)+4 9 (r)) _ u r \hi( r \i 

2b'(r)f(r) + + 8b(r)f(r)g'(r) + ^ fe - 8&(r)</'(r) 

- 8b(r)g'(r) 2 - 6(r)$'(r) 2 - 8e 2/{r) ~ 4s(r) = , (B.26) 
b"{r) + b\r) (/'(r) + 4#'(r)) = , (B.27) 
/' (r)/i'(r) + h(r) (2/"(r) + 2f'{r) 2 + 8g"{r) + 8g'{r) 2 + $'(r) 2 ) + Ag'{r)h'{r) + h"{r) = , 

(B.28) 

b'(r)f (r) - 4e 2 /M- 4 ^) , , 4 ^ Jf _, , „, , 2 



6(r) 



+ /" (r) + 4/'(r)</(r) + f (r) 2 = . (B.29) 



In the absence of any back-reaction (i.e. setting Nf = 0) we will recover the usual AdS- 
Schwarzschild xT 1,1 background given by 

r 4 

h {r) = — , b(r) = 1 - -f , $ = , (B.30) 



e / = e 9 = r . (B.31) 



In terms of the redefined radial coordinate p, defined by 



e -f- 4 °dr = dp e f +4 ^ = 4~ (B.32) 

dr dp 



25 These equations of motion are obtained by taking various linear combinations of the Einstein equations and 
the Klein-Gordon equation. 
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the equations of motion take the following simpler form 
b"(p) = , 

$" ( p) + _L &( p )b'(p) = 3( ^^ : ^ (p)+2/(p)+63(p) > 



9 " {p) + = bjpj 

h(p)h"(p) - //(p) 2 + ^ti(p)b'(p) + ^ = 

with an additional constraint equation, 
h"(p) 



_ 2e 4 /(p)+ 4 9(p) + 6e 2/(p)+6 9 (p) _ 3(2N f )g* e $( p ) +2 /(p)+6 g ( P ) 

47T 



16L 8 



$'(p) 2 - 24^' (P) 2 - 16<?W(P) + 2/" (p) + %"(p) = 



(B.33) 
(B.34) 

(B.35) 
(B.36) 
(B.37) 

(B.38) 



Note that the equation for b(p) completely decouples from the rest and is insensitive to the back- 
reaction. We can solve (B.33) and (B.37) analytically, yielding 

b(p) = b (l-£) , (B.39) 

h(p) = ±— j— cosh [ mp) feoj • (B.40) 

The constant pu denotes the location of the event-horizon. Also, bo, ho and h\ are three undeter- 
mined constants. 



Appendix C. The perturbative solution 

Recall that the equation of motion for b(p) completely decouples from the rest of the equations 
and can be solved exactly. We can find perturbative solutions for $(p), f(p), h(p) and g(p) by 
employing the following expansion 

$(p) = $o + e$i(p) , (C.l) 
f(p) = fo( P ) + ef 1 (p) , (C.2) 
9(P) = fo(p) + e 9l (p) , (C.3) 
h(p) = h (p) + ehip) , (C.4) 

where $ , h (p) and fo(p) represent the KW-background solution and it should be noted that 
e /o _ e 9o f or KW solution, given by 

ho(p(r)) = ^ , $o = 0, e kM) = r . (C.5) 
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Here we define 



6 *AT 3 ( XN i 



The redefined radial coordinate p is given by 



and h Q (p), f (p) are given by 



/ (p) = ~ln(-4p) 



Vp) 



-4R 4 p . 



Using equation (3.23), b(p) is given by 



b(p) 



1 



PH 



(C.6) 
(C.7) 

(C.8) 
(C.9) 

(CIO) 



We have set bo — 1 since this parameter just scales the time-direction. Therefore, at the first order 
in e, we obtain the following equations 



*i(p) 



1 



Kp) 



<&'i(pWp) 



1 



16p 2 b(p) ' 

fi(p) + T^-Mp) b '(p) + ^(p)) 

b{p) p 2 b{p) 



fx(p) + 7 9l (p) 



16 



p(pH - p)K(p) + (p - 2p H )h[(p) -h^p) = . 



Equations (|C.11[ C.14|) can be solved analytically, yielding 

P 



$i(p) =ci - j^ln 



c 2 ln 1- 



^i(p) 



c 3 



P# 



Ph 

+ 2(c 3 + 2c 4 ) + c 4 (2 



P# 



In 1 



PH 



(C.11) 
(C.12) 

(C.13) 
(C.14) 

(C.15) 
(C.16) 



Now imposing the boundary condition $i(p = p*) = will set C\ = 0, where p* represents an UV 
cutoff. The constants C2, c 3 and c 4 are still undetermined. The constraint equation (B.38), to the 
first order in e, becomes: 



2(/r(p)+4^'(p)) + 



K(p) 

h {p) 



16^(p)(/((p)+4^(p))=0. 



(C.17) 
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Using equation (C.16), we get 



2(f^p) + Ag'l(p)) + -(f[(p) + Ag[(p)) 



C4 



4L 4 (p-PhYph 



The last equation can be solved analytically to give: 

c 5 , ^(l-^)ln(l 



h(p) + ^gi(p) = c & + — i + 



p 

pn 



Eqs. (C.12 C.13) can be rewritten as 

flip) 



flip) + 

g'lip) + 



PH 



p- p H P 2 PH ~ P 3 



PH 



P~ Ph P 2 Ph - P 3 V 4 
Combining these two equations, we get 

1 , „, , , p H 



ApL 4 



(fi(p)+9i(p)) 



-hip) + \gi{p) - li 



if'M + ±9iip)) + 



P~ PH 



U'M + ^g'M) 



P 2 PH - P 3 



2(/ 1 (p) + 4< 7l (p)) 



1 

16 



Now using equation (C.19), we get the following condition: 

l(c A - 2c 5 ) . 



1 

32 



L 4 p H 



and flip) can be written as: 

(c 4 - 2c 5 ) 



flip) 



1 

32 



Up 



H 



\ C4 ( X 



p 

2p H 



In ( 1 - 

PH 



ApL A 



c 5 

pL* 



With that we have the following equation for g±(p): 

ig"{p) I 4g ' l(yP ^ I 39l i p ) pH - pH 



C4PH^Og (l 



P 

PH 



P~ Ph P 2 (p~ Ph) 32p 2 (p - p H ) Ap z L A (p - p H ) 

c 4 , c 4 log(l-- 



+ 



4p2L4 (p - PH ) 8p*L* (p - p H ) 
CbPh c 5 



(C.18) 



(C.19) 

(C.20) 
(C.21) 

(C.22) 
(C.23) 



(C.24) 



(C.25) 



fL* (p - pn) 2p*L* (p - p H ) 
Boundary conditions: 

Before solving for gi(p), we can determine some of the constants by demanding regularity of 
the Ricci scalar on the horizon. Using the equations of motion, the Ricci scalar S is given by 



p ._ p syMN 



bjp) 
hip) 1 / 2 ' 



-2(/(p)+4 9 (p)) 



itb{p) 
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,*(p)+2/(p)+6 S (p) 



(C.26) 



Now the condition that S is regular at the horizon tells us: 



c 4 = , c 2 = . 



(C.27) 



Note that the condition C4 = could have also been obtained from demanding regularity of fi(p) 
at the horizon. We can also choose c\ = 0. After redefining the constant C5 ( L ? H — C5), we have 



b(p) = 
$i(p) 



1 - 



PH 

16 VP* 



h 1 (p) = c 3 2 



/i(p) = ^(-16c 5 + l) + ^-% 1 (p) 



with (?i(p) obeying 



V/(p) + ^+ 3 " 1(PW 



Ph 



c 5 Ph 



+ 



CuPh 



P-Ph p 2 (p-p H ) 32p 2 (p-p H ) p 3 (p-p H ) 2p 2 (p-p H ) 
The general solution of this equation is given by 



(C.28) 
(C.29) 
(C.30) 
(C.31) 

(C.32) 



/ x 1 C 5 p H C 5 



0^^(1,1; 3;^) ^(2^(l 



g 



Ph 



Ph 3/2 \fplPn 
where, -E'(x) and iiT(x) are the Elliptic functions. We need to impose cj = for g\(p) to be regular 



(C.33) 



at the horizon. Therefore, using equations (C.31BC.33), we get 



9i(p) 



c 8 (lE ( 1 - - pP "< 

1 C 5PH C 5 8 V V *W 

"96 5p 10 ^/p7^ 



Ph 



f ( n 1 , C 5PH c 5 

/l(p)= -96 + ^p--10 + 



4c 8v /pK (l - ^) 8c 8v /p^ (l 



g 



v/P 



(C.34) 



(C.35) 



Now using equation (C.31) and the fact that dr = exp(/(p) + 4g(p))dp, to the first order in e we 
have 



r(p) 



and rn is given by 



V2(-p)^ 
1 



1 fe, 32- 2 + ^p-) 



V2(-p H y/ 4 



l + e 



1 3c 5 
32 10 



(C.36) 



(C.37) 
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Equation (C.36) can be inverted easily to get p(r) 



p(r) 



1 

'4^4 



1 

32 



C5 
2 



+ 0(e 2 ) . 



(C.38) 



Finally, we have the first order solutions in terms of the r coordinate: 



b(r) : 
$(r) 
h(r) 

e f(r) 



1 - 



H 



1 - 36C 5 

5 



e 
4 

L 4 



In 

+ e 
1 + e 
1 + e 



2c 5 + 



4L 4 c 5 



5r H 



+ c 3 2 - 



H 



2Cs 

5 

2C5 

5 



24 
1 

48 



+ 4c 8 



K 
K 



H 



H 



+ 2c 8 



H 



E II 



H 



(C.39) 
(C.40) 
(C.41) 
(C.42) 
(C.43) 



After absorbing the constant factor (l — |ecs) of 6(r) in all the other functions and then redefining 
the x % coordinates, the background solution can be given by 



b(r) 

$(r) 

h(r) 

e f(r) 



r 



e 
4 

L 4 



g 
r 



In 

■ + e 
1 + e 
1 + e 



- ) + + c 3 



2 
5 ( 



5r 



1 

"24 
1 

"48 



4c 8 



- c 8 



K 1 - 



K [ 1 — 



-8c 8 
+ 2c 8 



£ 1- 



-E 1 



(C.44) 
(C.45) 
(C.46) 
(C.47) 
(C.48) 



The solution presented in ( C.44 )-( C.48) is the most general solution of the equations of motion, 
which contains three undetermined constants denoted by C3, C5 and c 8 respectively. Below we will 
fix two of these constants and in the main text we will consider the background which eventually 
will have only one independent constant. 

In general, the solution for h(r) takes a form which is deformed compared to the warp factor 
which corresponds to pure AdS. Imposing the condition that h(r) ~ L 4 /r 4 we get 



c 3 



2L 4 c 5 
5r% 



(C.49) 
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Alternatively we could choose both C3 and C5 to be zero to begin with. This choice of constants, 
in the dual field theory, corresponds to setting a dimension 8 operator to zero. Thus we have fixed 
all constants of integration except c 8 . 

Now, it is important to note that the constant c§ can be a function of r#. In the limit r — > r*, 



e 2g(r) g 2/(r) 



r 2 + 0(1) (ec 8 ) 



2 ' 



(C.50) 



where r* 3> r#. But far away from the horizon, the metric should be independent of r#. Therefore, 
c$ is given by 

c 8 = r>" 2 . (C.51) 
The constant m -2 does not depend on r# and it has the dimension of length -2 . Finally we have 



b(r) 

$(r) 

h(r) 



e m = r 



1 



' H 

r 4 



e. 1 r 
4 ln 



1 + 
1 + e 

1 + e 



1 

"24 
1 

48 



4m' 



K 1 



m 



-2 



H 



K 1 



- 8m- 2 r 2 £ { 1 - if 



2m -2 r 2 £ 1 



yy 



(C.52) 

(C.53) 

(C.54) 
(C.55) 

(C.56) 



Appendix D. Fluctuations of the classical profile 

Let us denote the embedding of the probe brane by 



1" = + 



u a 



(d.i; 



where Xq denotes the classical profile of the probe, and denotes the transverse fluctuations. 
The indices /i, v = 0, ... 9; the indices i, j, k = 8, 9 and finally we will make use of the indices 
a,b,c = ... 7. Clearly, the indices fi, v etc. are used for the ten dimensional background, the 
indices i, j, k are used to represent the transverse fluctuations (which in our coordinate system are 
denoted by 89\ and 5<p\ respectively), and finally the indices a, b, c are used to represent the world 
volume of the probe brane. 



The background metric Gn V can be expanded under (D.I) as 



G, u = G% + d p G%E^ + -d p d a G%m^^ + ... . 



(D.2) 
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Here G Q denotes the background metric evaluated at the classical profile of the probe. 



Now, the induced metric (up to quadratic order in the transverse fluctuations) on 
D-brane can be obtained from the following formula 

9ab = G^d a X»d h X» 

12 
= 9ab + 9ab + 9ab > 

where 

9ab — G®d a XQ d b Xl 5^5 b , = G° d a X(}d b XQ , 



and 



gib = d p Gld a x^d b x^ + g% [daXgd&'q + wfl , 



g 2 ab = G%,dJ?d&>6?q + -d a X^d b X^d p d a G%^8^ 
+ d p G%^5f[d a X{;d b V5» + d b X»d a V5»] . 

Now the action functional for the probe brane is given by 

Strobe = -T 7 ^ d 8 £e®y/-E ab - J C 8 ^j , 
where £ represents the world volume coordinates, C 8 is an 8-form potential and 

E ab = El b + E l ab + E 2 ah , 

with 

J?0 _ E l _ 1 , F rp2 _ 2 

^ab — 9ab > ^ab ~ 9ab + 5 -^afe — i/ab > 



where F ab is the vector fluctuation. 
Now we use the following fact 

det(-E ab y/ 2 = det(-E° ab Y 

where 



/2 



1 + -TrM - -TrM 2 + - (TrM) 2 
2 4 8 V ; 



+ . . 



M c a = 9* El + g$El . 
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Now it is straightforward to check that at the quadratic order, we get the following contributions 



-TrM 



- (TrM) 2 

— -Tr (M) 2 
4 



n ab n 2 

% Sab 



1 n cb n a ' b> F TP n^n^b'n 1 n l 

■^9o % f ha >t h < c - -g g g a , h g Vi 



| [9o 9o - 9o 9o J 9b'c^ba' 



(D.12) 
(D.13) 

(D.14) 



Note that the above formulae are completely general for any background metric and any probe 
brane profile. So, in general, the vector and the scalar fluctuations will be coupled through the 
last term in (D.14). On the other hand, the Wess-Zumino term gives the following contribution 
at the quadratic order 



-P[C 4 ] A F A F 



(D.15) 



where -P[C 4 ] denotes the pull-back of the background C 4 potential and F is the vector fluctuation. 
Now the background metric we have considered is of the following form, 



ds 2 = h- x ' 2 {-bdt 2 + dx 2 ) + h 1 /' 

,2/ 



dr 2 e 29 



{dOi+ sin 2 94$) 



i=l,2 



6 " 2 

+ — [dip + cos 9\d(j)i + cos 2 d(j)2) 
9 



(D.16) 



The induced metric for the profile given by 6\ = it/2 and <\>\ = const is obtained to be 

(Jjl p 2g 



e 2f 

+ — (dip + COS 6*2^02 ) 

9 



Using the above informations and the expressions in (D.5)-(D.6) it can be shown that 

9* = (dJMt + dtSfaS*) = . 



(D.17) 



(D.18) 



The last equality follows from = cos9® = cos (tt/2) = 0. And 

9 2 ab = Gfc^a^W + G'^WlWl 

+ defile (dMA 2 + dtfifaS?) + de^l^e, (dj^f + (D.19) 



48 



Further simplifications follow: using the contributions listed in (D.12)-(D.14), it is also straight 



forward to check that for the induced metric in (D.17) 

1 



9oV 



(D.20) 



Hence the vector fluctuations decouple completely from the scalar ones. We will not worry about 
the vector fluctuations henceforth. 

Finally, the Lagrangian for the scalar fluctuation is obtained to be 



C 



~E° ab 



\gtG Q ei6l d a 8e 1 d h 5Bi + ^fGj^atyi Wi 



56t 



(D.21) 



The equations of motion for 89\ and 8(pi resulting from the Lagrangian in (|D.21|) are given by 

a, 



e^V^g^Gl^SOt 



(D.22) 



and 



e*V^E° (d^G^gt* + d^GP^ftf*) 56, 



(D.23) 



It is clear from (D.22) and (D.23) that in general the scalar fluctuations are coupled; specifically, 



they are coupled via the fluctuation modes along ip and 02 directions. However, using the expres- 
sions for the background and the induced metric, it can be shown that the coupled equations take 
the simpler form 



a, 



e^V^gfGl^h] - e*V=^fy$0i = . 



(D.24) 
(D.25) 



For the above equations the method of separation of variables cannot be applied. 
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